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/iBSTRACT 

The design of 4-'bar mechanisms for finitely separated 
positions of coupler guide has been considered. An analytical 
technique is formulated to synthesize for maximum of 3 or 5 
possible positions, depending whether the positions of the 
crank are coordinated or non-coordinated. When certain 
design and side constraints are imposed, an optimum solution 
is sought where inequality constraints are handled through 
Penalty Function* 



CHIPTER I 


INTRODUCTION 

Kinematic synthesis techniques have been developed . 
and refined to the point where^at least in the planar area, 
the kinematician has a fairly good idea of the type of motion 
that Can be obtained with a 4-bar linkage. Freudenstein jl--^ 
in his earlier papers, provided the background necessary to 
encourage the widespread interest in analytical kinematics 
that followed. With the invent of high speed computers the 
trend has been to replace conventional graphical techniques 
by sophisticated analytical methods where accuracy of design 
and economy of time are of prime importance. The present 
vjork has been carried on similar lines. All solutions have 
been obtained on IBM 7044. 

1.1. Link Guiding Problem - Definition ; 

Graphical solutions fs^of the 5 position problem with 
non-coo rdinated crank is quite well known. When the positions 
of the coupler are to be coordinated v/ith crank, the mazimum 
number of positions to be specified is reduced to 3, 

In all the above cases we get an exact solution of 
the 4-bar linkage. When certain other constraints such as 
maximum and minimum transmission angle , the approximate 
location of pivot points etc, are imposed on the problem, 
we may not in general, arrive at any exact solution, even 
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for the above cases. What \Te can do at most is to find the 
best solution •which will give a mechanism tha,t satisfies all 
the constraints and at the same time gives the minimum error 
for the positions when the crank is rotated through desired 
angle s . 

In the present dissertation, IJ-positions of the 
coupler and corresponding rotations of the crank are speci- 
fied'. In this case no closed form solution can be obtained. 
The best solution is arrived at by making a mathematical model 
and then seeking an optimal choice. For the given positions 
an analytical solution is sought for any 3 out of N-positions 
and this becomes the starting solution for iterations. In 
this method, the mathematical programming problem with 
constraints is transformed into a sequence of unconstrained 
minimization problems. 

1.2, Review of the Earlier Work : 

The concept of optimization has been applied to a 
n-umber of problems in kinematic synthesis in recent years. 

The work of Kiss j6j involves a generalized approach in 
terms of basic and characteristic motion criteria. Dynamic 
programming techniques are applied to solve the function 
generation problem. To solve f'unction generation problem by 
4-bar linkages, Garrett and Hall [?] randomly generated a 
large number of linkages and this was followed by the choice 
of the best solution on the basis of minimum root mean square 
error. Timkofsjalso attacked the function generation problem 



3 


in the same way. Fox and Wlllmert [sj have tackled the curve- 

I 

tracing problem in an elegant way. The problem of synthesizing 
a 4-bar linkage is presented as mathematical programming 
problem vjhero inequality constraints are handled through 
Penalty Function. Fox and Moore 10 have solved the same 
problem, the difforcncc being that they have assumed the error 
function as a 'continuous^ function while Willmert took it 


as a 'precision* point problem. The approach outlined by 
Han |\lj holds considerable promise for the application of 
dynamic programming to a wide range of design problems in 


kinematic synthesis. Gustavson |l2J applies weighting factors 
to three necessary design criteria to sort a large number of 


solutions, for point synthesis in terms of Burmester Curves 


for 4-finitely separated positions. Another unique approach 
for point synthesis in terms of damped least square iteration 
is given by Lewis and Gyory l3 .• Eschenbach and Tesan 
have solved the problem for 4-gene raliz;ed coplanar positions. 
The constraints are obtained from pairs of corresponding points 
on the generalized Burmester Curves,- These linkages are sorted 
by means of zones and given design criteria based on weigh- 
ting factors applied to various necessary design requirements. 
The present dissertation is an attempt to give an 



unified approach to point tracing, function generation and 
link guiding problems of 4-bar linkages. The approach is 
believed to be unique in the sense that this is the first 
attempt of this kind and most general in the sense that It seeks 
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the host design, for any given criterion, for, not only the 
well-known 3,4 or 5 position problems, but also for any finite 
number of given positions. 3, 4 or 5 point problems in 
point tracing, function generation or linlc guiding are but the 
particular cases of the approach presented. 

In Cha.ptor II analytical approach to the problem, for 
limited lumber of positions, us3ng complex numbers has been 
discussed. It has also been shown, that how the maximum 
number of 5 positions for non-coo rdinated crank reduces to 3, 
when crank is required to occupy the assigned positions. 

In Chapter III mathematical model of the problem in 
terms of objective function and design constraints is developed. 
Chapter IV deals with the solution algorithm* Few numerical 
examples have been Included in Chapter V. Discussions and 
Conclusions have been given in Chapter VI, 



CtlAFTER II 


lIIALYTICii approach TO THE PROBLEFi 

The method of complex numhers is applied for the 
kinematic synthesis of d-har link guiding problem. It is 
shown that upto 5 positions, closed form solution can be 
obtained for non-coordinated crank problem. In this case the 
designer has no control over the design. In 4 positions 
designer has the control over the design by varying arbitrary 
parameters. The method has been programmed for automatic 
ditigal computation on IBM 7044, and is included in Appendix. 

2.1. Coordinate System and Notations : 

In Figure 1 X-Y, is the coordinate system with ’0’ 
as its origin and P as the coupler point with as its 

coordinates. The linkage as shown consists of link 1, 2, 3, 4 
and 5 having lengths a, b, c, d and 1 respectively and angular 
positions 5 ©5 3} 05 H-"* 

2.2. Position Equations : 

Loop 1 consists of 1-inks 1, 2, 3 and 4 and position 
equation can be written as 

, 3© ^ 33 3 *=^ . ^3 4 * (9 i") 

be +ce' =ae +dem 

Loop 2 consists of links 6, 2, 5 and 7 and position 


equation can be written as 
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Since the motion of each link is to he expressed, in terms of 
its rotation from axi unspecified starting position, it is 
convenient to write 


0 = 0^ + t 

,S> = a + u 
f 1 

V + u 

and0 = 0^ + V 

Writing equations (2,l) and (2.2) as 


he e + c e ’ 


1 e3^= ae^^+de^^lq.^^ 


f 0 ^^ + h e 


33t ^ 1 ^3^1 


= r e^^^ 


(2.3) 

(2.4) 


Wow if we define vectors as 

.3 ^ 


1 


■■6 


= a e' 

= b e^®l 
= 

= d aJ'*! 

= 1 '+^1 
.3 A 


= f e- 


R 


r e 


3>^ 


Hence the position equations can he written as 


Zg e 


3t 


4“ Zo G 


- Z 4 = Z. 


(2.5) 


= R 


7 j. 7 ^ 3 t . „ 

Zg + Z2 0 + Z5 & 


( 2 . 6 ) 
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2.3. Uso of Position Equation for Synthesis : 

In Figure 3 various positions of guide are given. To 
specify the location of a guide a.t any position, its orienta- 
tion and the coordinates of any point on it, say coupler point, 
must he known. Since the guide is fioted on the coupler, its 
orientation in the i"^^ position from original position would 
he Uj^, as that for coupler. Hence the problem is defined hy 
specifying coordinates of coupler point (X^, Y^) and u^. . 

In the form of equation 

“^0 ” ^i ( i = l, 2 ,,.,,jn) 

Therefore, the problem for (n+l) positions could bo defined as 
Zg + Z 3 e^^^ - (i = 0,l,...,,n) (2.7) 

Zg + Zg + Zg e^^i = (i = 0 , 1 , ...,,n) ( 2 . 8 ) 

Solution of equations (2*7) and (2,8) gives the 
required synthesis. 

2.4, Determination of iiaximum Number of Coupler Positions for 
Coordinated Cpank Problem ; 

In this Case the guide is required to occupy various 
specified positions corresponding to given angular positions 
of input crank. 

♦ 

Known quantities for (n+l) positions a,re 

, R^ for (i = 0 , 1 , . ... ..,n) = n+l 

tjL for (i = 1 , 2 , ,n) = n 

i 
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and unMiown quantities are 

for (i = l 52 ,,.,,,n) = n 

Zi 5 Zg, Z 3 , Z^, Zg and Zg 

Since each vector requires two real quantities for its speci- 
fication, its magnitude and direction, tho six vectors give 
12 unknowns. Hence total number of unknowns are 12 + n* 

Each - complex equation gives two real equations, hence 
total number of available real equations, referring to equa- 
tions (2,7) and (2.8), are 4(n+l) . Hence the maximum 
possible positions for which the problem could be solved is 

4n + 4 = i2 + n 
"max = 2 

Therefore maximum number of positions for this case would be 
3 ( i = 0,1,2) , 

2.5. Determination of Maximum Number of Coupler Positions for 
Non-coo rdinated Crank Problem ; 

In this case number of known quantities are 

Uj; for (i = 1 , 2 , ,n) = n 

and number of unknown quantities are 


^i 

for 

(i^l jS ^ • 1 


= n 


for 

( i"”l ^2 5 , 

• • • jH) 

= n 
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Six vectors- giving 12 unknowns. Hence 

the maxinium possihle positions for which the problem could 
bo solved is 

4^ + 4 = 12 + 2n 

^ax ” ^ 

Therefore maximum number of positions for this case would be 
5 Ci — Ojlj...» ^4) ■ 

Now onwards, 3 position problem shall refer to 
coordinated crank problem and 4 or 5 position problem to 
non- coordinated crank problem. 

2.6. Synthesis of Three Position Problem ; 

As explained in section 2.3 wo can write equations 
similar to equations (2.7) and (2.8) for 3 positions 

Loop 1 ; 


^2 

"^3 

’ \ 

= 

= 



(2.9) 

^2 


+ Z3 

3^1 

e 

- 4 

e'^l: 

= Z 

T_ 

(2.10) 

^2 


■^23 

,3U2 

- 4 

e3V2 , 

= z 

1 

(2.11) 

Loop 2 : 












jwo 

0 



(2.12) 

^6 

^2 



3^1 

e 

“ ^1 

•r-D 

CD 

(2.13) 



3*2 


3^2 


awg 

(2.14) 

^6 

+ Z2 

e 

■"25 

e 

= ^2 

e 
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It is clear from above equations tha-t there are 14 
unknowns ,2 for angles v^, Vg and 12 for 6 vectors, but we have 

only l2 real equations at our disposal. Therefore in order to 

get a nontrivial solution wo can keep v^ and V 2 arbitrary. 

First wa solve equations (2.i2 - 2,14) for Zg, Zg 
and Zg and then solving (2.9) - (2,11), we find Z^, Zg and Z^. 

A Careful examination of the above equations reveals that no 
solution of the problem exists if (t^ = u^ and t2 = ^ 2 ^ and/ 
or (u^ = v^ and Ug = Vg) . 

By subtracting equation (2.12) from equation (2,l3) 
and (2.14), vje get two equations with no Zg term 

an Z + b Z^ = c (2.15) 

ag Zg + bg Zg = Og (2.16) 


where 


a^^ = 1 - e 


bi = 1 - e 




3^i 


(i=l, 2 ) 


c^ = ro e ^ - r^ e 


(2.17) 

( 2 . 18 ) 

(2.19) 


^2. 


Multiplying equation (2.15) by (-* ~) and adding equation 
( 2 . 16 ) : 


where 



_ -K ^2 

t-S- 



( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 
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Knowing Z from equation (2.20), we can find Z from equation 

O Ci 

(2.15) and Zg from equation (2,i2), 

Writing equations (2.9)-(2.1l) as 


- Zs = Z2 


2 ■t'Z e^^“2 ^ - o ^ 1 


1 


Zg e 


2 ^ 2 . 3 to 

Z^ + Z^ 0 - Z 3 e = Zg e ^ 


(2.23) 

(2.24) 

(2.25) 


These equations are alike (2.i2)-(2,i4) , only difference being 


a. - 1 - e ^ 
1 


bi = c 


= p3^i - 


1 


(i = 1 , 2 ) 




3t. 


(2.26) 

(2.27) 

(2.28) 


and Zg, Zg and Zg are replaced by Z^, 2 ^ and Z 3 < 


2,7. Synthesis of Four Position Problem : 

We Can write position equations for this problem as 


Loop 1 : 


Z, + Z^ - Z^ = Z„ 
1 4 3 2 


2+2 e^""! - Z e^"^! = 


1 4 


Zo e 


3t 


1 


\ + 2^ ” 23 e^^ = Zg 


3t, 


2 + z . 2 = Zg 


_ 


ot. 


(2.29) ' 

(2.30) 

(2.31) 

(2.32) 
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Loop 2 ; 


^6 **" '^2 ^5 ~ ^0 ^ 


^6 ■" ^2 ■*■ \ 

Zg + Z„ 0^*2 ^ ^3^2 

u ^ 5 

Z + 7 e^^^ + Z 

6 2 5 


3W. 


r e 
1 


^2 


^3 " 


3^2 


3^3 


(2.33) 

(2.34) 

(2.35) 

(2.36) 


It is clear from above equations that there a,re 18 unknoT/jns, 

6 for angles t^, tgj t^, v^, Vg and v^ and 12 for 6 voctors, 

but we have only 16 real equations at our disposal. Hence to 

get nontrivial solution we keep 2 quantities as arbitrary. 

First wo solve equations (2.33)-.(2.36) for % 

and Zg, In this set we have 9 unknowns, 3 for t^, tgj t^ and 

six for 3 vectors Zg, and Z^ but only 8 real equations. 

Hence we choose t^ arbitrary and express tg and tg in terms 

of t . 

1 

Proceeding as in section 2,6 V70 have \ 


^1 ^2 ^5 

= "l 

(2.37) 

^2 ^2 ^2 ^5 


(2.38) 

^•3 ^2 ^ ^3 ^5 

zz ro- 

(2.39) 

- 1 ^3^1 

a. - 1 “ e 

1 


(2.40) 

b^ = 1 - e^^i 

( i = 1, .... ,3 ) 

(2.41) 

c = r . 

1 ^0 ® 

' ^i ® 

(2.42) 
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And, 


^2' ^6 “ 




(2.43) 


lAThe re 


bi ' = b 




(i = 2,3) 


= c-i- 

^ 1 ai 


ajU 

_ h . t t 

% 

Using equations (2.43) - (2,44) we can find characteristic 
equation. 


(2.44) 

(2.45) 

(2.46 


^ 2 ’ ^ 3 ' 


(2.47) 


c^a*, b/^a-i Coa-1 "^oa-i 

or (c ~^) (b = (c - (2.48) 

1 a^ 1 ao 1 a^ i- 2 


1 a3 


or (Cgbg - C2b3)a + a2^a2 ^ ^ 1 ^ 3 ^ ^ 2 ^l” ^1^2^”^ 

(2.49) 


Defining d^, dg and d3 as 


°3^2 " ^2^3^ s-l 


'^2 “ ^3^1^ 


d-n - (^ 2^1 " a. 


1^2^ 




(2.50) 

(2.51) 

(2.52) 

(2.53) 


Writing equation (2,49) as 


- dg e 


i'^2 = rl. e^*3 


(2.54) 
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Taking complex conjugate of botli sides 


-T -r "3^2 - -jtq 

- d^ e = do e 


(2.55) 


Multiplying equations (2.54) and (2.55) and rearranging terms 


^ 4,^4 ^2*^2 *" ‘^-3 '^3** ^2^4 ^ - "^2^4® 


“3t2 


0 


(2.56) 


Defining d^ as 


<^5 = 'i4'i4 + d2'^2 “ 
Equation (2.56) can be written as 


(2,57) 


dod^ e ^ - dc- e ^ ^ + d^d^ = 0 

4^ 41: O 4 


(2.58) 


This is a quadratic equation in e and its two roots are 
given as 


3^2 _ % - ^*^5 ~ ^^ 2 


3t2 

1/2 


2d2 d^ 


(2.59) 


3 "^2 3 "^3 

Knowing e , e can be found from equation (2.54) as 


Vs _ 14 . 1 ^ 2 .! 




(2.60) 


itp jts 

Knowing e and e we can solve for Zg, Zg and Zg using 
equations (2.43), (2.37) and (2.33). 

The set of equations (2.29) - (2,32) is alike 
(2,33) - (2,36). There are 9 unknoims, 3 for v^, Vg and Vg 
and 6 for 3 vectors Z^, Z3 and Z^. Hence v^ is chosen 


arbitrary 



3.7 


Whe re , 

1 

= 1 - e ^ (2.61) 

1 (i = 1,....,3) (2.62) 




(2.63) 


and ve can proceed in the same \>:ay for solving Z^, and Zg. 

A carefu.1 examination of the above equations reveals 


that no solution exists if (t^ = u^, tg = u^ and t^ - u^) 
and/or (u^ = v^, u^ = v^ and u^ = Vg) . 



CHAPTER III 


M;jrHEiviAIIC.iL F0Ei4ULiiII0I\T 


In this section, the mathematical model of the 4-bar 
linkage is formulated. The term ’ma,thematical model* means 
the complete set of equations, and ineqnalities, which is 
necessary for the statement of the design problem as a mathe- 
matical programming problem. This set of functions includes 
all equations needed to describe the geometry of the 4-bar 
linkage, the motion of its coupler guide, the objective 
function, and the design constraints. 


3.1. Objective Function : - . 

The object of the synthesis problem being considered 
here is to desi.gn a 4-bar liniage, where coupler guide follows 
as closely as possible some assigned positions, and whose 
crank rotations are as close as possible to a set of desired 
values. The desired positions are given in terms of coordi- 
nates of the coupler (x. , jj-) point P and the orientation of 

JL 

the guide . 

Since guide is fixed to the coupler 


ts.n . = - ■‘I. 

* 1 vl H 


.C. f. 


1 a 

The desired rotations of the crank are 


C i"~l ^2 5 • « ♦ • ^ n—l) C3#l) 


AY. = T- 

■ 1 •* 1 


+1 


C i— 1,2 j * * * . ,n— l} (3.2) 


I 
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Consider the d-har linkage shown in figure 4. 
Coordinates of coupler point are functions of linkage parameters 
(X ^5 , . . . and of input angle Xj;i,(='^i’ ). The 9 

quantitie s (X^, . , . , . ,Xg pCg, . , . , ^X^q) are talcen as design 
Variables in the synthesis problem. We define a design 
vector D as 

Tx 

1 i 


D = 


X, 


“•8 


^10 


(3.3) 


Any set of values for 9 design variables, D is 
called a 'design', even if design is absurd. The coordinates 
of coupler point may bo written as 


2^= X (D, f 25***‘5 


-> n 


7o~ J (^1 ^o?****? 


(3 ,4) 


(3.5) 


(3,6) 


.p ^ ^ .^n' 

and orientation of the guide may be written as 
A’] = -VKd, 5 Xg,...,, j^) 

An error function may be formulated as the sum of 

three functions, B^(D, . . . . , T ^) ,’ ’ * * * ’ '^ n^ 

So(D, T which define the coordinate matching 

error, crank angle matching error and guide orientation 


matching error. 
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S'tG:.rtlng with the initial design oi a 4-bar linlcaga , 
the cranh is rotated through given angles, each t.lme solving 
the linkage for coordinates of coupler point and orientation 
of the guide in terms of the design vector D and and thus 
formulating two error functions instead of the three, as 
mentioned above. Thus, tho total error function is 

E ( B , ■^c +"^1 2 Eq (D, 

(3.7) 

where, (j~ and (T” are weighting factors chosen to gi’n the 
12 

desired relative importance of two types of errors. 

flow including in design vector B, we have 


s (d)='T; e (b) +cr' e^(b) 

1 c 2 0 


(3.8) 


The 'best' design is the design ’..diich minimizes the 
tota.1 error, E(D), and s^-tisfies all the design constraints. 

What now remains in tho development of the objective 
function is to determine the coordinates of the coupler point 
and guide orientation as function of design vector B. 
Befarring to figure 4 

O 1 /o 

(3.9) 


p 2 T_/2 

B = (X^ + “ 2X7 Cos (X^)) 


1 


T 2 


T X9 Sin (X7) 

A = sin'^ 

D2+X42/X32 

( gJiT" ) 


) 


<1 


- 7\ 


(3.10) 

(3.11) 

(3.12) 


X 


3x 


= X^ Cos (Xg) + X^ Cos + X^) - X9 Cos (X.7+ (3.13) 


^1 


"4 


8 ' 


"8- 
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Sin (Xg) + X^ Sin,<§.+ Xg) - Xg Sin (X^+ Xg) (3.14) 
-1 Xq 

The coordinator of coupler point can be written as 


Xp = Xg Cos (X^q) + Xg Cos Xg Cos (Xg +^j) 
Yp = Xg Sin (X^g) + Xg Sin (X^^ Xg Sin (Xg +'>| ) 
Now the total error function is 


E(D) = -J 


1 


n 

Y" 

i=i 




.-1 f 


2 X ' 


(Xp - x) + (y^ - y) +(Y 


i=l - 

,2 


i 


(3.15) 

(3.16) 



(3.17) 


3,2 Design Constraints s 

The design constraints, a satisfaction of which 
distinguishes acceptable designs from unacceptable designs 
Can be grouped into two categories* /uiy constraint _which 
restricts the range of design variables for reasons ether than 
the iirect consideration of performance of the design is 
Called a side constraint. L constraint which is derived from 
the explicit consideration of the performance of the design 
is called a behaviour constraint. 

3.2,1 ’ Realizability Constraint : 


We may say that all lihlc lengths should be greater 


than a given minimum values 
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h- 

X . 
min 

> 

0 

(3.18) 

''2 ■ 

X 

min 

> 

0 

(3.19) 

"^3 - 

^min 

> 

0 

(3,20) 

^4 - 

X . 

min 

> 0 

(3.2l) 

Also requiring that 

^7 

should he greater 

than zero and less 

than 360 degrees give 



X 7 $ 

1 0 



(3,22) 

2V - 

Xy ^ 

0 


(3.23) 


3,2,2 Crank Rocker Constraint : 

Req-uiring that the 4-'bar linkage he a crank- rocker 
mechanism adds another set of five constraints. The Grashoff 
Ine qual it y|l9j states that if the sum of the lengths of longest 
and shortest link is less than the sum of other two links, and 
if the shortest link is the input link, then the 4-har linkage 
will he a crank rocker mechanism. This gives 


:i0 

(3,24) 

X3 - Xg ':^o 

(3.26) 

X 4 - Xg ^0 

(3,26) 

X3 +X 4 -Xj_ -X 2 ;^o- 

(3.27) 

(Xj^- - ttg - X 4 )®;^o 

(3.28) 


3,2,3 Compactness Constraint : 

Since there are usually limitations on the span in 
which the mechanism is to he used, a set of constraints is 
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used to limit the length of each link to less than a prescrihed 
mazimum length. It is not necessary to place size limit on 
link because the length of X^ is limited by constraints 
placed on the location of the pivot points O^^and Og. 

■^maZ ~ ^2 (3,29) 

^maz *“ ^3 ^ (3,30) 

3.2,4 Pivot location Constraints i 

In some mechanism design problems, it may be 
necessary to limit the pivot points 0^ and Og (figure 5) to 
lie in prescribed regions of the XY plane. For simplicity 
a circular area is chosen here, however it could be taken 
to be of any definable portion of the plane. 

The limitation of the location of 0^ gives 

fx, - X.! R, (3,32) 

' ii 6{ kn ii. 

or in terms of design variables 

^A ^6 ° (3.33) 

and the limitation of the location of Og 

But 

Xg- Xg- X3_ ® = (Xg- Xg- X^) ■ (Xg- Xg- xp 

= Xg^ + * Xg® + SXg Xg Cos ^8^ 

-2XgXgCos(4-Xg-Xgo)-2XBXgCos (^Xg-Xg) 


(3.35) 
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Therefore the inequality (3.42) may be written as 
%^“^B^~^6^"^l"**2K3XgCos(X^Q-X3)+2X3XgCos( 4 

+2XgX4os( /v Xg-Xg) >0 (3.36) 

where R^, Rg, X^^j Xg 5 >^ Xg can be chosen by the 

designer, 

3.2.5 Force Moment Constraints : 

There is limitation of the torque and forces in the 
links as a result of the force applied on the coupler point. 
The axial reaction force in link 4, Kg in link 2 
and the torque in the input link can be expressed as 
Appendix b. 



(3 ,4l) 
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! 


gW^SinCO^-X 



2 


X5 ¥j_Sin(ni4-XQ.>Q.) 
• X3Sin(Xgj_ - h ^ 


:>0 

(3.42) 


3.2,6 Transmission Lngle Constraints : 

Transmission angle should he as close to '/7/2 as 
possible for better transmission characteristics. Maximum 
and minimum transmission angles occur when crank is in line 
of base link. 


. = Cos“^ 

i min 


• ,(. = Cos 

: max- 


-1 


3X3 X4 

XgS + X42 - (X4 + Xg)^ 
*“ 2X3 X4 


(3.43) 

(3.44) 


One common practice is that the deviation of andM*-.^ 

• maz * mm 

should be less than 30®. 


and 


' ; ^ I tj ^ 

6 *' ■ 'max 



; -TT 



Ka. f/O 

mih 


(3.45) 

(3.46) 




CHAPTER IV 


SOLUTION ALGORITHI4 

The method employed in the solution of constrained 
minimization prohlem is discussed in this Chapter. The 
constrained problem is first converted into an unconstrained 
minimization problem and then solved as a sequence of uncons- 
trained problems. Each of the unconstrained minimization 
problems in the sequence is solved by the iterative descent 
method in which the gradient vector is approximated by 
diffe re nee s. 

4.1 Conversion to Unconstrained Minimization : 

The geperal mathematical programming problem is to 

determine' a vector X that minimizes tile function f( X ) 

' >1. 

subject to the constraints gj_CX) //O , ( i = 1,2, . . ^m) 

Fiacco and McCormick! 20,2]^ have developed an alogo- 
rithm for transforming this mathematical programming problem 
with constraints into a sequence of unconstrained minimiza- 
tion problems. The procedure is based on the minimization of 
a new function 

m 

p (X,r) = f(X) + r ^ (4*1) 

i=l ^ 

over a strictly monotonic decreasing sequence of r-values^;^;^* 
Choose r£> 0 and a point Xq which lies strictly 
within the constraint set, i.e* gi(XQ )^05 (i = l,2,,«,,,m). 
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It is expected that a minimum of p (X,r^) exists inside the 
constraint set 5 since on the boundary of the constraint set 
some of gj_(X) must equal to zero, and p(X,r^) The 

minimum point of p(X,r-,) depends on the value of r. , and is 
denoted by X(r^). The point XCrq^) lies inside the constraint 
set. By reducing the value of r, the influence of the summa- 
tion term, which penalizes closeness to the constraint boun- 
daries, is reduced, and, in minimizing p(X,r) more emphasis 
is placed upon decreasing f(X) If the minimization process is 
repeated for a decreasing sequence of r values, r£^ r^'^ 

r_- \ 0, each minimum point XCr.^) can move closer to the 
boundary of the constraint set if it is profitable in terms 
of reducing f(X). This method is particularly attractive for 
problems, such as the 4-bar linkage problem, which have 
nonlinear constraints, since it approaches the solution from 
points inside the constraint set, and thus, avoids the diffi- 
cult movement along the nonlinear boundary of the set, which' 
is required in other methods. 

Fiacco and McCormick have shown that if 

(i) the interior of the constraint set is nonempty, 

(ii) the functions f(X) and gj_(X) , (i=l,....,m) are twice 
continuously differentiable, 

(iii)the set of points in the constraint set for which f(X).<VQ ' 
is bounded for every finite vq, and • 

(iv)the function f(X) is bounded below for X in the constraint 
set,- 
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then, the optimal solution to the unconstrained problem 
a.pp roaches a local minimum of the constrained problem as the 
value of r approaches zero. If, in addition 

(v) f(T) and gj_(X), (i=l,2, . . . . ,m) are convex functions, and 

(vi) p(X,r.) is strictly convex in the interior of the cons- 
traint set for every r 0, 

then, the optimal solution to the unconstrained problem- 
approaches the absolute minimim of the- constrained problem as 
r approaches zero. 

In order to apply the Fiacco McCormick algoritlim 
to the 4-bar linkage problem, the constraints must be expressed 
as penalty functions. The penalty function for the constraints 
in section (3.2i) to (3. 25), aro simply the reciprocal of the 
constraint functions, since these are all discrete constraints* 

4.2 Choice of Initial r i 

To start the algorithm a decision has to be taken 
about the initial value of r and factor c, the factor by which 
the value of r shall be decreased. The algorithm, imposes no 
restriction on the value of r anc C except that r^ 0 and C )> 1« 
With largo r- values the number of computations required to 
reach the constrained minimum increases tremendously while 
with small value of r, the function is highly distorted and 
requires initial solution very close to actual minima* 
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4,3 Method of Unconstrained Minimization : 

jin unconstrained m-inimization of the function 
p(Djr) must he carried out over the design vector, Dj for 
each value of r used in the Fiacco -McCormick algorithim. 

The method of unconstrained minimization used in this 
dissertation is a powerful quadratically convergent 
iterative descent method, orginaliy given by Davidon[22j, 
and subsequently modified by Fletcher and Powell(23, 24j'. 

The method is based on the properties of the 
quadratic function 

0 (X) = 0 + X^b + ^ X^ (4.2) 

where 0 ,,,. is ascalar, (t- , ...... ) , b is a 

constant n-vector, and n is a ssnnmetric, positive defi- 
nite matrix. The gradient of the function is 

g = b + iiX ' (4.3) 

If we consider the quadratic form (4,2), then, 

0 

given the matrix = tkX-lSx • » calculate the 

1 3 

displacement between the point X and the minimum Xq as 

Tq- X = - A“^g • (4.4) 

In this method is not ovaluated directly. Instead a 
matrix H is used which may initially (H^) be chosen to be 
any positive definite symmetric matrix (identify matrix I). 
The method is descent method which proceeds towards the 
minimum of the function by a sequence of linear minimi- 
zations. The directions along which the function is 
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minimized are chosen so that if the function being minimized 
is quadratic function (4.2), its minimim will he found 
after n linear minimization, 

i. typical iteration in the minimization of a 
function f(X) proceeds as follows ; 


1, Given the starting point Xj_ and the gradient of the 
function evaluated at gj_ = f (X^) , the direction 
along which the (i+l) minimum will be sought is calculated 
from 


ii = - [Hfj ii 

Where H is positive definite matrix. 


(4.5) 


2, FindT'j so that f(X^ +'^. d^) is a minimum along the 
line Xj_ -*4" dj^ 

s. set 

4, Calculate the new gradient vector -17 f (X t 4 .;]_) 


and set y^= g^+-L- gj_ 

5, Calculate the new H matrix by 


- _T 
di di 


- T 


L^i+i'i = i tT - 


d. Xi 


r.®-il ^i ^i L ^^il 

- T - 

Xi ^i Xi 


6, Begin the next iteration from step 


(4.6) 


The modification of the matrix (step 5) is 
done to take into account local characteristics of the 
function in order to avoid the zig-zag characteristics 
common to many other minimization techniques, such as the 



34 


lELstliod of steopest descent. It can be shown that the 
matrix H^l approximates the inverse of the matrix of 
second partial derivatives of f(x) evaluated at the 
mininnjjn because as the minimum is approached the second 
order terms predominate in Taylor's expansion. It can 
also be siiown that if the function being minimized is the 
quadratic function (4.2), then 


In case minimum along a line is not found, a 
new [h]) matrix is not calculated in step 5. The new 
direction is calculated from equation (4.5), using the 
new gradient, gj_+ 2_5 and the old matrix. This always 
helps to ensure that f H ] is always positive definite, and 
thus the directions dj_ will always be directions of descent* 
It is found that is no longer positive definite, then 
Hj_ is reset to the identity matrix, '“l], and the iteration 


is continued. 

The unconstrained minimization procedure is 
terminated when the components of both direction dj_ and 
the actual stepT~. d. become less than a prescribed level. 


4.4 Gradient t 

In order to calculate the directions used in 
Fletchor and Powell's method, it is necessary to calculate 
the gradient vector of the function f(X) during the course 
of each iteration. Since analytic expressions for the 
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partial derivatives of p(X r) are quite cumtiersome, the 
gradient vector is approximated by differences. 


Let 

and 


f 

f* 


r f (At j . . . 


'll’ 

”*’^r***‘* hn 

= f* 



Then, by Taylor’s series expansion and neglecting 
order terms 


(4.7) 

(4.8) 
highe r 


_ f*" f 

Z - ^ -5 ^ t 


(4.9) 


S3 

A simple scheme ife follovedto calculate as., v. can 

be approximated as .whesCfeS, is a small positive 

quantity. Each time when f is changed to f*^ all constraints 

are checked to ensure that none is violated. If not.Af ^ 

' 'I J 

should be replaced by ^^ 4 * Sven if all constraints are 

^ 2 

not satisfied should be replaced by Xj_. Thus all 

components of the gx^ldient arc calculated. If any^^ 
turns out to be zero, corresponding should be 

replaced by a small positive quantity • 


4.5 Linear Minimization ; 

In the Fletcher and Powell procedure it is 
necessary to find the value of -A = such that fCX^^+'A^d^) 
is the minimum of f(X) along the line X = X^ +a( d^. The 
method used for this minimization along a line is Golden 
Section SearchfBd'j. Since the Fletcher- Powell technique 
is based on the sissumption that f(X) can be approximated 
by a quadratic function in the vicinity of the minimum, 
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"bho liiiGsp ni.lni.nii . '3 s.'t ion scIigiiig is TDcisod. on "th.© niini!iiiz&~ 
tion of a quadratic function. 

If and d^ar© considsrsd fixsdj f(X) can b© 
considered as a function of only, 

f(Z) = f(X^ +Tdi) = f(;r) (4,10) 

Then,'.,}' ^ is the value ofl" for which fC^r") is minimum, 
LetX^i^ he the value of'X beyond which constraints are 
violated. 

First we start with some positive value of >(, and 
check whether replacing X^ by X^+^ d^^ all constraints 
are satisfied or not. Thus, we determine th© bounds ofc< 
aS'X„.,„ for which constraints are satisfied and'-^^^ _ for 
which constraints are not satisfied. 



f(c>^) as After determining Golden Section 

search can again b© used to find the minimum of fiok ) as 
fC'-Af) between’!7<,= 0 and'-.K.= ^^^ 4 * 




Reduced Interval Reduced rlnterval Reduced Interval 

iof uncertainiy i of uncertainly i uncertainly 
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G-olden Section Search proceeds as follows i 
Considering the -animodal function f(X) in the 


Interval a ^ h evaluated at two intermediate ordinates 
as shown in figure 8. In the case shown in figure (8 .a), 
the extremum lies in the interval In the case 

shown in figure (81) the extremum lies in the interval 
a <^X ^ X 2 and in X s^X 2 for case figure (8c), 

In this procedure and | are so defined that 
if the region a^d'X^'X^ is eliminated ^ the ordinate X 2 
is in the proper location; that it is the first ordinate 
of tho new intervaljA-X = b - X^ i»e. located att^AX 
from X } "the new left end of interval ef uncertainty} 
figure 10. Similarly, if the region Xg /^X was 
eliminated by X^ being greater than Xg, then X^ would be 
in the proper location to be the second ordinate in the 
interval^ X^ = (X 2 - a) i.o. located at >1 from a. 

This readiness for either alternative demands 
a symmetry in the location of X^ 'and Xg which gives 
:!= 0.381966011 
] = 0.618033989. 



chapter V 


TTa^iERIC.J. EXi’jjyiPLSS 

In this Chapter few numerical examples have 
heen solved, 

5,1 Three -Position Synthesis : 

10,9, 8.4, 6,0 

y 11.7, 13.0, 5,2 

pi 



40,0 deg. 


-19,0 deg 

*2 

80,0 deg. 

U2 

-20,0 deg 

Arbitrary 

Pararaeter 



"'l 

5.0 deg. 



^2 

25,0 deg. 



Solution 





32.8 



^2 

4.7 



CO 

15 ,2 



X4 

24,1 



^5 

5,5 




14,1 



T 

60,4 deg. 



4 X 

6,2 deg. 



4x5- ^3 

97,0 deg. 



4x6 

18.6 deg. 




* 
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5,2 Four-Position Synthesis : 
Desired Motion 


y . 11.7, 12.8, 11.9, 9.7 

u^ -12.0, -19.0, -20,0 deg. 

Arbitrary Parameters 
t^ 40,0 deg, 

7,0 dog. 

Solution 


4 

4 

X4 

4 

^6 


^ ^5” 

iae 

'^2 

^3 


32.0 

5.4 

2.6 

34.7 

19.3 
32.5 

205.0 deg. 

147 .4 dog. 
“18.4 deg. 

94.7 deg, 
7,2 deg, 
3.1, deg. 
38.1 deg. 

11.3 deg. 
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5.3 Five-Position Synthesis : 
Desired Motion 


^pi 

10.9, 9,6, 7.3, 5.3, 5,6 


11.7, 12.8, 11.9, 8.3, 6.0 

^i 

-12,0, “19.0, -19,0, “7.0 deg 

Solution 

X 

12.0 

1 


4,0 


9,5 

X. 

8 . 0 , 

4 


4,0. 

5 

a 

X 

8.5 

6 

j 

31.0 deg. 


11,0 deg. 

3 X 5 -XX 3 

29,0 deg. 


36.0 deg. 
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5,4 SevGU-Position Synthesis - Optimization : 
riote .s All angles arc in degrees. 


Desired 

Motion 


X , 

pi 

10.9, 

9.6, 7.3, 5.6, 5.3, 5.3, 5.6 

Xpi 

11.7, 

12.8, 11.9, 9.7, 8.3, 7.1, 6.0 

^i 

-12.0 

, -7.0, -1.0, 1.0, 6,0, 5.0 


40.0, 

40.0, 40,0, 40.0, 40.0, 40.0 


Constraints 

1. Crank- rocker mechanism 

2. Xg, Xg>/0 

3. X = 20 *0 5 X — 1»0 

max nin ri 

4. = 10. 0, 10.0, 10,0, 10,0, 10.0, 10.0, 10.0 

5. ©. = 10.0, 10,0, 10,0, 10. 0, 10.9 5 10.0, 10,0 

6. = 4.0, X^^ = 8.0 and = 30.0 

Rg = 4.0, Xg = 20.0 and^Xg = 25,0 

7. Tg = 50.0, K 2 = 20.0 and = 30.0 

Solution 

1. Value of C = 4.0 
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Design Variables 

Initial 

Final 


Solution 

Solution 


12.0 

11.9492 

^2 

4.0 

3 .9783 

^3 

9,5 

9.5015 

^4 

8.0 

7.9859 

^5 

4.0 

4.0320 

^6 

8.5 

8.5015 

r 

V 

31.0 

24.8317 

4 

11.0 

13.1213 

-'-4~Xg *• 

29.0 

19.1214 


36.0 

34.8319 

r 

10.0 ^ 

.0015 

E(D) 

8.2061 

1.0664 

p(D,r) 

123,5217 

1,0664 



CHi;jPTER VI 


CONCLUSIONS iUID DISCUSSIONS 

In the analytic approach to the problem it is found 
that 4-posttion problem will find wide spread usage, as the 
designer has the control oyer arbitrary parameters and can 
obtain variety of mechanisms. Though the 5~position problem 
will give moro accurate results for a given N-position 
problem, but, it is complicated one. In this case wo have 
to solve a polynomial of 8th order which is really a 
cumbersome job. 

In optimized solution it is known that conditions 
(i)“(iv) (Section 4*l) are satisfied for the 4-bar linkage 
problem. However, since 'there is no simple way to test the 
convexity conditions, (v) and (vi), for this nonlinear 
problem, it is not known whether the optimal solution to 
tho unconstrained minimization problem Is the optimal 
solution to the constrained minimization prohlom or only 
a local minimum of the constrained problem. It is expected 
that the method will tend to converge to a minimum which 
has the greatest Influence in the vicinity of the starting 
point. Therefore, the final design for a given 4-bar linkage 
problem solved by tho method presented here is dependent 
tho choiCG of initial design. 


upon 
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The main problem which arises in using the 
Penalty Function Approach is the choice of initial value 
of Fiacco and. IicCorird.clc have presented two methods 

for choosing the initial value of ’r', both of which 
depend on some knowledge of the characteristics of the 
gradient vector. In using the finite difference 
formulae to compute the gradient at each step, these 
necessary factors are not available. In the present 
work Initial 'r' is taken to be ’lO’ which gives pretty 
good results. Two methods of finding initial 'r’ are 
given in Appendix 

The gradient calculation which has been done 
in the present work is .not very accurate. A better 
scheme is given by G.¥. Stewart, III for estimating 

the gradient vector by difference. 

Other areas of further study include investi- 
gation of the present problem wdth design constraints 
such as excluded region ’for coupler motion, dynamic 
characteristics and material property* 


G.W. Stewart, III, "A Modification of Davidon's 
Minimization Method to Accept Difference 
Approximations of Derivatives , oournal of 
the Association for Computing Machinery, Vol.l4, 
No.l, January 1967, pp 72-83, 
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APmiDIZ A 

SS^TRESIS OF FIV15 POSITION PROBLEM 
le Can write position equation for this case 


clS 

Loop 1 : 

^4 ■ ^3 ” ^2 

Z + 

14 o 

^ 4 ® “■ ^ 3 © 


Z + z .e ^ = 

1 ^ O' 


jV3 

+ Z^e ^'^4 - 236^’^ = 
Loop 2 : 


Zge 


Zge 


Zge 


Zge 




jt2 




0t4 




3^0 




Ze + + Zge^^^l = 

Zg + Zp0^^2 + 256 = 


Zg + Zoe'^''^3 + z„e'^^ 


r,e^-l 


rge. 


3^2 


^5®’ 


rge 


3^3 


3^4 


Zg + ZpG + Zge 




It is clear that there are 20 unknowns, 
t^,**»*, "^4 e*nh xT"^,*#*** 


(1.1) 

(^.2) 

(A.3) 

(A.4) 

CA.5) 

■ 1 ^. 

(A*6) 

(A.7) 

(A.8) 

(A.9) 

(A.lO) 

8 for angles 


,v^ and 12 for 6 vectors. There 
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TG al G quations at our disposal* Hghcg a nontri- 
vial solution Can "bo oltained uriquoly without assigning any 
parameter arhitrary. 


where 


Procooding similar to section 2*7, we have 

ai Z2 + Zg = (i = (A.ll) 



i' (i — 1,»»»» }4) 


°i = 




7^0 


- r^ e 


DWi 1 


j 


(A. 12) 


(A,l3) 


(A. 14) 


From set of equations (A.ll) we get 

!)!_ = Cj_* (i = 2, ,4) 


who re 



(i = 2,... ..,4) 


(A. 15) 


U.l6) 


Eouations in (A. 15) give the following two independent 
equations s 


and 



(A.17) 


(A. IB) 
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Expanding equation (^*.17) in terms of and using 


following identities 

di - 031^ - cgbg (A. 19) 

dg = c^b3 - 03!^ (ii.20) 

d3 = Cgb^ - c^bg (A.21) 

d4 = -(d^ + <^2 (A.22) 

VJg get 

d^ o^H + dg + d4 = - d3 (.^^23) 

Multiplying by its complex conjugate 

d^e^^l + dgS’^'^S + d4 = - d3e^''3 (^1.24) 

we got 


^ '*^1^ + d d e^'^l+d d + 

12 1 ^ 41 14 

d^dg + d4d2G^^^ = ^3*^3 (A.25) 

Defining 

p = d4d2 (A, 26) 


q - d^d^ 


s - ^^^2 


r = dg/S 


(A.27) 
C A«28 ) 
(A. 29) 
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Equation (A, 25) can be written as 


;i(t2-ti) - _Kt2-ti) .X _ 

pe +pe +q + q e ^ 1 + s e 

- -;it2 

+ s 0 + 2r = 0 

Sixiiilarly equation (ii«iS) can be written as 

D(t2-tq) _ -jCtp-t.) - -DtT 

P' e -^ + p' e ^ 1 + q* e ^ 


Dt2 


(A.30) 


+ s'e ‘^ + s*G ‘^ + 2r»=0' (A.3i) 

whore following identities have been used 

d ^ “ ^^.'^2 "" ^ 2*^4 Ciik * 32 ) 

5-2’ - ~ (a. 33) 

dg ^ — ^2^]_ ^q^2 (n«34) 

d4’ = (d^‘ + dg' + dgO '(A.35) 


p’, q’, r’, s' are similar top,q.,r,s as in 
(A. 26) - (A,29), where d^ is replaced by 

dq' . 

Recalling that e = Cos x + i Sin x, we obtain 
real parts of ocuation (ix.30) and (A.3i) , respectively 
as fellows, 

Px ^■^2’’'^1^ " Py ^'^2"'^1^'^ ^ ^1 “ V 

+ s^ Cos tg - Sy Sin t 2 + = 0 (A.36) 

p^' Cos (tg-tq)- Py' Sin (tg-tq) + Cos t - qy»Sin tq 


+ s^’ Cos tg - Sy' Sin tg + r^-’ =0 (A.37) 
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These two uquations are real and independent, transcen- 
dental in the Wo unknowns, t^ and t 2 > To transform these 
into polynomial equations, we utilize the following trigo- 
noLiGtric identities. 


Cos 9 


Sin 9 


1-tan^ i 
1+tan^ i 

„ 9 

Stan 2 

l+tan^ I 


(A.38) 


(k.39) 


To suhstitute those identities in equation (A.36), we 
first simplify the notation hy letting 


Tu = tf.ui ^ (3 = 3,4) 


(^.40) 


Thus , e qua i'ion (A. 36) hecomos 

I 1-^1 ^ l-Tg^ 4T^Tg ■' 

i + T^2 ^^^^2 + 





+ = 0 (A.41) 

Equation (A,37) will give similar result? the difference, 
however will be the presence of primes (’). 
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Hq'^ation (ii«4i) gives 

+ C^) * Tg(,.^T^2 ^ Bgl^ ^ Cg) 


+ (AgT^ + BgT^ + C3) = 0 
Similarly equation (A. 37) will yield 

^ qCigT^® + Bgl^ + Cg) 


(a.42) 


+ (A_^ + C^) ^ 0 (A.4S) 


wiie re 


^'*1 - Px - ^ - Sx (il-44) 

— — 2(py + q^.) (il*45) 

= (^ - + r,, - (ii.46) 

"'^3 " ""x “ Px “ qx 

(ii.»48) 

*^3 = Px \ 

■^5 “ “^Py “ ^y^ (i^.50) 

Bg = 4p^ (A,5l) 

Cg = -2(py+Sy) "■ (A.52) 


Exrrcssions, for A, B and C wi-th the subscripts 2, 4, 6 
are in terms of p^’, c^', s^’ etc., and are similar to 
expressions shoim previously. 

'pplying Sylvester's dyalitic method of elimi- 


nation, we liiultirly both equations (A*42) and (A.43) by 

2 3 

Tg and consider T 2 and T2 as different variables. 

This yields a system of four linear equations with unknowns 


Tgs "^2^ ^2^ 

of unknown Tq. 


(e.g. X, Y, Z) with coe fficients in terms 

f iii.iiKmmr I 

f MiTi'i? .. i^isiAay* 
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In order to got a nontrivial solution, we must 

have 


0 


AglihBsTi^S 

AoT.^+BoT^+C, 

0 

2 

A 2 T ^ +B pT ^+C 2 

2 

AgT^ +B0T2^+Cg 



Vl^+Bsh'^S 

aIoT *1 "hb oT -i Q 

0 j_ 0 ^ 

0 


VlhBgTi+Cg 

A4Ti^+B4Ti+C4 

0 


(A, 53) 


This elimiriejit is the corupat ibility equation, which 
must be satisfied by T^. 

This determinsiit set equal to zero and expanded 
wj.ll yield an eighth degree polynomial equation in T^* 


KgT^®+ iqi P.gT/+ HsT 3_®+ Rj^T3_+ Rq^O 


Where E.'s are given as follows : 
D 

r=8 = h-'i-E" 


R? = 

Rg = D^ng+DjR^+DgSg-aF^Fg-Fg® 

Eg = Di 34+D4E3+D3E2^D2S3'*2 FiF 4“2F2F3 

= D^E5+D5S^+D2S4+ D4E2+D333-2F^F5“2F2F4-^3' 


(A. 54) 

(A.55) 
(A. 56) 
(A, 57) 
(A.58) 
(A.59) 
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=■' °2’^6+I'5'’2+I’3S4+D4l3-2FgF5-2y3F4 

(A. GO) 

Rg = D33j.+D4E4+D5Rg-2F3F5-F42 



(A.62) 

Eq = DgEg-Fg® 

(A. 63) 

¥ht.'ro D's, 3';-: and F's are .'■jiven as follows t 


D ~ ^ • 

(A. 64) 


(...65) 

^3 = •‘*4^5'^*'‘5’^4'*'^4®5’''‘'^3'^6“^‘6*^3"®3^6 

(A. 66 ) 

= ®4^5‘^®5^4“B3C6"^6^3 

(:*.67) 

’^S = ^4^5-^3C6 

(A. 68 ) 

■p A /- n • 

(A.69) 

Eg = A^Eg+ngB^-i^Bg-xigEg 

(A,70) 

33 = 3^Cg+.lgC^+B3^Bg-;.2C5'-35C2- BgBg 

(A.71) 

3^ = L 2^0 qC ^*^620 g — B gC 2 

(ix.72) 

5 16^0 

(A. 73) 

?! = iYl^~_ig ..2 

(:..74) 

Fg = -*-1®4''*-^’^4®1~-‘‘3®2‘"'^^2^3 

(A.75) 
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F 


3 


^4 ”■ “® 2*^3 


^5 ^ 1 ' 4 "^ 2^3 


T' 

3^2 


(...76) 

(i..77) 

(...78) 


This equation yields 8 roots which are complex pairs and/ 
or real roots. For each real root of we can easiljr find 


t 


1- 


ts - 2 tan ^(T^) 


(1.79) 

To obtain tg , ue return to couationG (1,42) and (;.,43). If 
those have common real root(s) for T^, then 
tg = 2 tan ^ (Tg) 

To find tg and 14 vie use 


(A.80) 


tg- Log ; 


1^1 ^ ^ ^2 ^ ^4 


. 1 , 

^4= ^ Log 


1 j-t 


e 


1 ^"^2 
•^+ dp^ e + d4' 


-d. 


t 


(1.81) 


(1.82) 


_j 


2g, Zg and Lg can he found from (...15), (;..ll) and (1,6). 

Similarly Z4, and Z^ can he found as Zg, Zg and 

Zg, only difference heing 

DU, 


1 


Ui = 1 - e 

■u 

hj_ = e -1 


Cf - Zg - Zg e 


V 

f 

{ 


(i = 1,4) 


(n.83) 

(1.84) 

(1.85) 
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APP3NDIZ B 

FORCE kOMSNT COlCSTPu-JH'ES 

Free Body dir..gram 03 ? the co’apler link is shown 
in Figure G, 

The sum of the moments about point must be 

ze TO . 

. 

0 ■, Xg Sin (Xg+'^'l) 

-W Sin 0 . Xg Cos (Xg+ n) 


- jK^ Cos i-rr - ? - Xg) Xg, Sin V|1 


- IK 4 Sin ( Tf - Xg) Xg . Cos , T]! 

(B.l) 

i^i 

0 = Xg. V. Sin (Xg+rj* 0 )+X 3 .K^.Siii(V| _£ - Jg) 

(B.2) 

Xg.w^.sin x_gi - ep 

■41- Xg.sLi (Xg^+S^ -np 

(B.3) 

IL 

K 4 ' - 2 _‘ 4 V 0 

(B.4) 

The sum of forces in the direction of must be 

equal to 

ze ro . 


W •Rq Kp ^ < 

^ _ ^2 ^ Kg Kg 

(B.5) 

0 = Kg+W Cos(X 7 +X 3 ~©)+ K 4 Cos(Xy- 6 ) 

(B.6) 

^ 21 = - ^41- 'h '®1' ^7i‘ ^8i> 

(B.7) 

n 

Kg? - 2ir k ^4 >.o 
i=l ^ 

^ ' If 

(B.8) 

f 
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The sum oi forces in the direction of vector ‘t' 
must ho equal to zero. 




I2 + ^Lt 

^2 (?-2) 

^2 


+ 




(B^9) 


^^2 

0 • Cos 

2 



- Xg + 0) + e^.Cos(A-, ^ ) 

(B.IO) 


The voctor t is defined such that X X t gives the iiioHient 

Tg* 


hl= ■ ^2 • ^4i- - ^71^ - ^2- '^1- 




n 

- xL t2 Q 

i=l 2i ^ 


(B.ll) 

(B.12) 



■# 


t 


^ 0 . 




# 


i 


t, 
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APPENDIX C 


C.'"OICS OF INITIAL VALUE OF 'r‘ 


Two iiietbods to choose the initial value of 
have been cuc-estod in the litorature. 


T 


( 1 ) 


j-^exo) 


HI 

where p(X^) = Z,.^ — p^-r 

0 g. (I^) 


1 ^ 

i3,l 

and represents the gradient of the function. 


(C.l) 


(C.2) 


CASS I : r^ / 0 

If r^ </ 0, minimization of f(XQ) alone, without 
considering the renalty term, is carried out. At every 
nex^ roint tlio r-value is checked by equation (C.l), If 
the value is positive, unconstrained minimization with 
obtained value of r is carried out. 

CASS II ; r^ = 0 

This means that the unconstrained minimum has 
been reached Xq = X , 

CASE III ; r^") 0 

If r^'^O, this is taken as the starting value 
for usual minimization. That value of r^ is selected 
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which ma].-oj objective function and Penalty term equal 

(C.3) 


£(X_(^)_ 


'1 


§ 1 ( 1 ) 

i=l 


111 

5 


.1 ]iis lias two weaknesses* Firstl^.^'j the starting 
point may ce too dose to some of the boundaries thereby 
mal ing such scDlection of r useless* Secondly j there may 
not bo a point in the feasiM.o domain which is not near a 
constraint boundary. If some of the constraints are near 
the boundary, i.e. penalty is large, the value of »r* obtained 
from (C.3) is very small. This difficultj’’ is eliminated if 
in the summation for penalty, these constraints which are 


quite close to the boundary arc substituted by some minimum 
value of constraint modified expression for r^ is 




f(X ) 
m 


(G.4) 


1 


[g^Xo), Vn 


Since constraints are normalized here, the value of 
can be arbj.trarily selected as 0,2. 
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ilOB MEG134 , T I MEOOa » PAGESO 1 / 5 N AME S , K » CHATURVED I 
SEXECUTE WATFGR , rewind 


c 


c 

c 

c 

c 

c 

c 

c 

c 


OPTIMIZATION PROGRAMME FOR FOUR-EAR LINK GUIDING PRCBLEM*^^ 

XP(I) corresponds to desired X-COORDINaTE of rOUPLER PO?Nf 

YP(I) corresponds TO DESIRED Y-COORdINaTE OF COUPLER POINT 

W(l) corresponds to force applied at coupler point at an angle THFTA(I) 
WITH the horizontal m t m.iMiij 

DETAII) gives successive angular increments in coupler positions 

DGAMA(I) g^ves succ-ssive angular increments in crank 

i''i gives Number of design variables in design vector 

NiN gives number OF CONSTRAINTS 

NP05 GIVES number OF POSITION SYNTHESIS 

FEASIELC REGION FOR THE PIVOT A IS DEFINED WITH RADIUS RA AROUND THE 
POINT FIXED BY RADIUS VECTOR XA WITH ANGULAR POSITION ALA 
X8 » RB, AND ALB ARE ANALOGOUS TO'XA* RA» AND ALA FOR PIVOT A 
X(I) TS THE STARTING SOLUTION OF THE PROBLEM 
COMMON/PI/PI 


COM'-'ON/MANE/N.' 

COMMON/iNEQLTVcONOO) tSATSFY 

C0M|''':0N/SUB1/DETA(6) ,DGAMA(6) ,XP{ 7) ,YG(7) 

C0MM0N/SUB2/XLMIN,XLMAX 
COMMON /SUB 3 / NPOs , N I N 
C O M M 0 iM / S U B 4 / s ' I G M A 

CCMMON/SUB5/W(7 ) ,THETA(7) , RA ,XA , ALA » RB ,XB , ALB , DESRT 2 , DESRK4 , DESRKZ 
SIGMA=1. 

DIMENSION Xf lO) 

PI=4,^(aTAN( 1. ) 

N = 10^ 

NIN = 26 
NPOS = 7 
NPOSl=NPOS-l 

^4T4 DETA/ — 12* ^ / o >6* j5« / 5DGAMA/4 d* >4q* ^4go )>3q*j2q**2q*/ 

^'ATA XP/l0*'?T9«S,7.3,5.6,ii.3,5.3,5.6/,YG/llo7,12.3,11.9,9.7,8.3, 

1 ^ o 1 » 6 , „ / 

Data X/l'-«!i4«t9«3»3»ji4»i>L»5)31«jll»»ii9»»36»/ 

DATA W/7«in, /, thETA/7*1o./ 

data XA»F;A,ALA»XB,RB,ALB,DESRT2,DESRKA,DESRK2/8. ,4,,3c. ,20. ,4., 


1.25 . ,30* ,2o* 9 3o« / 

ALA = ALA'^PI / 130 . 
ALB=ALn*Pl/ 180 « 

DO 28 1 = 1 , NPOS 

28 ThFTA( I )“THFTA( I )^-Pl/ 130 , 
DO 26 T=l,MPOSl 
DFTA( T )=DETA( I )*P 1 / 180 . 

26 DgAMAI I )=DGAMA( I )-*fPI/ 18 o. 

DO 27 1 = 7 , lo 
27 X( I)=X(T )^fPl/ 180 . 

Print so 

PRINT 100 ,(XP( T ) , 1 = 1 , NPOS) 

print 51 
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PRINT 100, (YG( I ) , I=1»NP0S) 

PRINT 52 

PRINT 100, (DGAMA( 1) ,I=l,NPOSl) 

PRINT 53 

print 100, (DFTA ( I ) , 1=1 ,NPOs1 ) 

XLriIN = l, 

XLMAX = 20 
Call .•iriMA(X} 

50 FORmaT(10X,-5;-XP( I )^<-) 

52 FOR-IAT ( lox ^^f-D-Af'iA ( I ) -> ) 

5 I FO RHAT ( ] ox 5 -'A YC- ( J ) ) 

53 format { 1 CX .*DrT A ( I ) * ) 

100 FORMATd CX»8F15^A) 

STOP 

End 

SIBFTC MINIMA 

Subroutine ■ (x i 

common/'.amf/n 

COMMON/M IN I /R 
COMMOI-'i/POVEL/Z 
COMMON /SERCH2/F UN VAL 
C 0 iM M 0 N / I N F Q L T / C 0 i V ( 3 0 ) , S A T S F Y 
Dimension x{1o),xi{10) 

R=io. 

281 1=0 

52 Call f(!nt(x,Z) 

IF( Id!F.0)G0 TO 53 

print ]'^0,FuNvAL 
Call Gutputcx) 

Print 70 , (con(J) , 0=1,235 

70 F0RMAT(21X,10F11oA-» 

CALL CHFCL (X) 

IF (SATSFY,>rQ ) oO TO 53 
PRINT 59 

59 FORMAT(//2nX,-“-lNlT9AL SSLUTlO.X dOfS M^T Ljf Vix CONSTRAINT 
STOP 

53 PRINT 600 
PRINT 5^'''' 

PRINT 650 ,R 
PPT.lT Si^'O 
Call'' i-C'.'FLL(x? 

1 = 1+1 

I F ( I c =0c 1 )C-0 TO 60 
DlFF=APS (Z-Zl ) 

IF(DIFF.LEo-''-BS(Z "i«E-Oo) )G0 TO 90 
80 Zl=7 
■- = R/4<, 


GO TO 52 

loo FORMAT! //10X,*T 

F 0 R r ' , A T ( 6 X , F- R , 9 X , ■<- 1 T R N , 


i I T I AE FUNCT I on value* ,E l5.6,l0X ,*SOLUT I ON 
5X ,*SFCCNDAHY* »6X ,M0 b0fcT I VE * , 


iTi.ON VfcT0R-«-/26X,-"-FUNCTICN-»,7X,*FUNCTI0'- 


SET*) 


VECTOR*) 

20X,*SOLU 
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c Value of alpha =c in linear mjn, means no progress in linear mjn. 

IF( ITN,EQ.^}G0 to 920 
IF (Z.Lt»Z1 ) go to 920 
C .THIS FOLLOVJS F { XA+ aLPhA^S ) = 

Z = Z1 
RETURN 

920 DO 92 I=i,N 

SIG‘-'A( I )=ALPl!'S( I ) 

ET=AMAW1 ( le sAE-'-S ( ALPI ) ) 

DO 76 L=1,N 

IF(ET^LS(L) „LF,.?PSL0N(L) ) GO TO 76 
GO TO 89 
76 CONTIMmE 
RETURN 

89 D0TSIG = D0TSIG+5IGMA ( n**2 
92 X( I )=X( I )+SIGFA{ I ) 

Call gradix^gd 
AR = 0 
BR = 0 
CR=0 


DR = 0 


DO 93 1=1, N 

ar=ar+gi n )«Gi( I ) 

BR=RR+s( 15*5(1) 

DR=DR+g1 ( I )*S( I ) 

.Y{n=Gl(I)-G{I) 

' CR=CR^Y( I )*S{ I 5 
93 G{I)=G1(I) 

DO 86 J=1,N. , 

ROC J)=R0(J) + (1./(aLPi«cR)-DR/CALPi*CR)**2 )*Y{ J)**2+2./CR*(''U)*G 
K J) ) 

IF{RO( J) cLE.O. 5 GO TO 71 
86 CONTINUE 


95 


98 

99 


A1=0. 

A2 = 0. 

DO 95 1 = 1, N 
A1=A1+SIGMA( I )*Y( I 5 
DO 95 J=1,N 
A2=A2+Y( I )*H( I ♦J)*Y( J) 
A1=1,/a1 


A2=1./a2 

DS 98 !V=1.N'^ _ / , 

Dd 98 j\l *N V , { . 

Hix I ,J)=W( I ,J)+S^IGMa\i )*^gKAi( J)*A1 
DO 98 K=1,N 

DO 98 M=1,N 

Hl( I , Jy=Hl( I ,J)-«-A2*H( I ,K)*Y{K)*Y(M)*H(M,J> 
DO 99 r=l,N I 

DO 99 J=1,N f 


HcI,J)=H1{I,J) 

itn=itn+i 


66 


PRlflT 5r>ro, ITm»Z ,F 
Call output (X) 

IF( TT:'^ara.J.'>)i-f.!) no jO 106 
Z1 = Z 

GO TO on 

106 PPir-'T lOT 

:c7 F0- ’AT( ///lox,''' 0f.xir.-vu;- f'.o OF ITFP '.T TL:.;5 

jnr'.r) FCP"' AT { 16y I ? , 2v , 2Fl 5 « 6 ) 

105 RFTUPP 
E p 

FlOFTn 

SUPPF.'UTTNF nRAD{X,G) 

CCiVC.i/MANE/.’-' 

COM-iOO / SUB 3 / NPOS ^ '1 1 1'' 

CuK'-Or /TRFCLT/r0f!(30) »5ATSFY 
CC; 'T)H/POVf L /2 
D I '-.el S T OG X ( 1 » X 1 ( 1 0 ) » f ( 1<t ) 

Fx = 7 

DO 100 J = l->N 
DO in 

10 xi(.i)=y(J) 

Dx=x ( T 1 
IF(rx*FQ,o») 0X=oOL 

5 XI ( I ) =X( I )+DX 
call CHECK(Xl5 
IF (SATFFV. FG. 0« } GO TO 20 
IF ( OX. LT»o,) GO TO 11 
DX=-0X 
GO TO 5 

11 Dx=-nx«„] 

GO TO F 

20 Call F"FT(v] ,f^'1) 

1^0 G( T )= ( rxl-- X 1 /OX 

retool 

EK'D 

GIEFTC F:P''x0*' 

FMOCTIO • "w;VXO'.(X'f 
CO;.'- C;VS'.i34/SK''.'. 

CGf ■•o:. /PI/PI 
CO • 0 /"AOE/.'' 

CO:. 0 ./IGF’-LT/C0;.;(30) TSFY 

CC. O'./SUp. i/O^T A ( 6 I , (0) tXP{ 7> .Yn{ 7) 

cc.- ,o- /cu:;^/ '’Cs,.:r.- 

C r . ' C- / A , '1 n / r T ft , r A *■.' a 
COF /SUB <=/ I 
U T ; ■ . r K, f j 0 f l X ( in) 

FijKXO" =n. 

Gaya --X ( 7 ) 

DO 1'-'' T=1.NP0^ 

CALL /'•GLF(X) 

GA'--A=G/':'-G- + X ( 3 ) 

ETA = FTft + X(‘'') 


PEOCF GUIT.!?-) 


Xa = '' ( 6 ) »CC5 f X ( !'> } ) + X ( 2 ) -^rOS ( 5 +X ( 5 ) -^COS ( ETA ) 
Ya = X( 6 ) ;;ST\ ( X « IC) )+X { 2 ) •'■ST:; ( G^Xat+X ( 5) *SIN{ ETA ) 
FUf'.XOi‘; = FUNXO!''i+ { Xa-X-'^ ( I ) ) YA-YG ( I ) ) **2 

GAMA = 6A(-,A-X ( ) 

EtA=FTA-X {*'’ ) 

IF<T,FO. 1)^0 T'; 90 

dftai=fta-etai 

IFT'^'ETaI ,Gr.2,'<^Pl )nETAl = DETAl-2.*Pi 
IF(GFTAl.Lf , (-2, ) x-Pf ) D’^TA1 = DETa1 + 2,»Pi 
FUi'''XO".' = PUN:- Of, +5 if.iA-;' ( DFTaI-OFT A ( I-l ) ) ^*2 
90 ETA1=ETA 

I F ( T . EO.NPO<; ) GO TO lOO 
GAMA = G aMA+OGA V ■' ( I j 
ICO CONTINUE 

return 

end 

SIBFTC FUNT 

Subroutine fumt(x,fx) 
com: O'-i/UINT /R 
C OKU, ON / SU B 3 / NPO S , N I N 
COM'VON / I NFQLT /con (30) , 5A TSFY 
COMT'O:^ /MANE/''-' 

C 0 M.VOM / SERCH2/FUNVAE 
Dimension xdo) 

FX=FuNxON (X ) 

fun\/al=fx 

CALL CONSTR{X) 

DO 10 1 = 1, ''■UN 
10 FX.=FX+R/C0N( 1) 

RETURN 

EAjD 

SIRFTC CHECK 

Subroutine chpck(X) 

COM’'-'ON/MANE/N 
COM' ,0N /SUB3/NP05 ,-N I N 
COI-.';CN/INFOLT/CON(30) ,SATSFY 
C 0 Mi'-.O i 4 / SU B 2 / X L I N , X L M A X 

C0M''’0N/SUBi/DETA(6} ,DGAMa(6) ,XP(7) ,YGC7) 

Dimension x(1o> 
satsfy=c 
Call cO'OSTt^(X) 
uG 10 I=l,NJN 
IF{CC i{ n «L>;^.0» )G0 to is 
10 continue 
ret< irn 

15 SAT"F/=1 
RETURN- 
E .A; D 

SIBFTC ANGLE 

subroutine ANGLE(X) 
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COMMON/ANr./FTA ? GA''';A 
CCMMOM/Pl/Pl 

C i M I ^ J 

0i,f;rrvSTO:''i xd-’) 

COfv. .O'-'./SUBS/v-, (7 ) ,THETA( '/ ) , RA , XA , f L A » R J , XiJ , AL B , DESKT2 , DESRK4 ,OESRK2 
C0i''"0 •./rONST/ T2,K.4,K2 
REAL GGjKA 

;<A=f' 

K2 = 0 
T2 = '^ 

D = X ( 1 ',-i: ;;2+X ( 2 ) ■:"«2-2 , '-^'X ( 1 ) ( 2 ) ’^COS ( C-AL-Ia ) 

D:= 5G.p_T f ) 

T r ( X ( 2 ) 1 R ( C AEiA ) / n ) 

Y L ['■; n = A Q c, I f T P ) 

CC ( ns.'-'+x { A ) ’-X ( 3 ) >'-5!-2 ) / { 2 , *D*X ( 4 } ) 

BET/ =A'-''CS(C'; «BET) 

D F L T A = F’ I - : r ■ T A - X L K Q 

X3X = X ( 1 ) " CLiS { X ( o ) ) +X { '+ ) *rOS ( DELTA+X ( b ) ) -X ( 2 ) *rOS ( GAi-lA + X ( 8 ) ) 

X3Y = X( 1 ) --GS 1 .dXi o ) )+X ( 4 )*Slr'.( DELTA+X ( 8) 5-X ( 2 ) *S I'A ( GA 'A + X ( 8 ) ) 

IF ( AnS(X3X) .LF. 1.E-3B )GO TC IqG 
Trv,p=AnS(X3Y/x3X) 

TE'--'-’=; t,aN(TE iP) 

IF ( (X3x.GT,0s, ) o AGP, (X3 y,LT,0,))TE -.P^-Ts^ ■•P+2 . *P I 
IF { (XA-X.LT.O. ) ,ANn, (XBY.CT.O, ) ) TT: .P=P I 
IF { (x3y.LT = n. 5 * ANo. ( <.3 y.LT ,0o > i T‘^‘'P = P I+TFf'P 
EtA=TEL'P 

i<.4 = .<.'.4 + {X{3 ) I ( ETA+X ( V )-TH'^TA { I ) ) ) / ( X ( 3 ) *5 I N ( X ( 8 )+DFLTA-ETA) 

1 ) 

T2 =T2-x ( 2 ) .:-K'V*3 In ( iiFLTA-GA. -X ( 2 ) :>v, ( I ) “ 5 1 ' ( THETA ( I ) -r-AR A'X ( 8 ) ) 

N2=|:.?-KAJicOT. (PELT A-CA' lA )-v, C ! ) i<COS ( TF’ET-'' ( I )-GA;-:a-X ( 8 ) ) 

RETURN 

100 IF (X3YoLT.f.« )rT'^=-Pl/2, 

IF (X3Yd:Eo-''e |-t,.-=xPt/2, 

A'.-, 4+( X ( 3 ! T )*SI ! ( ~T' +X( Y )-THPTA ( I ) ) ) / { X ( ' ) «-.S I N { X ( S ) +DFLT A-ETA ) 

1 ) 

T2 = T2-X { 2 )*A'- A CD'^LTA-^-.A-A )-X( 2 ) ( ll ’^-S lU ( T»^ETA ( I )-GA.‘iA-X( 3) ) 

.•.2 = 3 2- • '4 ■ C''5 ('■■'■LTa-G-''-:'.M-U{ T ) -^('C S ( T^FT a ( I )-GAMA-X ( 5) ) 

Rft: 

END 

TGPTC CO'XSTf' 

SUBPO'.iT! F C T.' ■■ TO ( X ) 

CO'T 0../i^;F .LT/Ccp (30) sSATSFY 
CO' 0 ./'n/T-'i 

c c. :: N / f3 u rj 2 / y L I , I i L . . A x 

CC . 0 /SUPS/NPC ..:C 

c. - . •Oi./OUdA/'; {'i ) •. THETA ( <' ) , R A > XA , A L A » RB , XB » ALB, DESRT2 » DFSRlC^, DESRt^2 

•■iFAL X2CK4 

Dli-NGTri'. 

COY C'vCOyST/ T2,K4,K2 
print, (Y( I ),I = -Ld.-) 

DO IC T=i,A 
cou( I )=x( n-xLd.-; 
lo CON( 1+4) -XLH.AV-X ( n 

CC0i(9)=X(3)+y {4)-y(l)-X(2) 

CON { 10 ) = < X ( 1 ) -X ( 2 ) ) "-^2- ( X ( 3 ) „X (4 ) ) ^*2 

DO 15 T=7,10 
CON( 1+4) =X ( I ) 

15 CON( 1+E ) =2CTPI-X( I) 

CON{ 19)=X( 1 )-y { 2) 

CON( ) = y ( 3>-X( 2) 

CON{21)=X{A)-y(2) 
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• CON { 22 ) = Rr- -2-X ( 6 ) -*2-X ( 1 ) *-'^2-XB*^^2-2 , ^:-X Q ) *X,( 6 )*cuS ( X ( 10 ) “X ( 8 1 ) 
1 + 2 o^^XB+X C 6 ) +rOS ( ALR-X ( 10 ) ) +2 e +XB*X { 1 ) +^0.S ( aLB-*X (8 ) ) 

COX' { 2 ■■ 1 =R a-Jh; 2-XA “ “B-X ( 6 ) *■" 2 + 2 ,+X A*X ( 6 } ^i-COS ( AL-A-X ( lO ) ) 

RBTU'' ' 

END ' 

-TC OUTPUT 

SijnROUTT'iF OUTPUT (X) 

CO 

0 I" ''~NstO+ X{1o) 
print !'■ ,(X(T) »T = l,N) 

0 FOP '.AT; ^bX-; S'-'l-'ob/////// ) 
aFTUP 
END 

FTC pFTA 

o'.JRROliTlNE pet, A {X,S, ALPHA) 

COi'i.'.’ON/CODf/ICO 

COl ON/SERC hI/F' ;lo»Fr.;AX»E!'1ID 

COri' Of'-'/f-iAfiE/N 

CC'-ooN/SUB/NTN 

C0i'';^0 '•■/INECLT/CCo (3o ) ,SATSFY 

Du'+'MstOM y ( lO) ,.S{ l''^) .y( 1 0> 

EtA = 2 « 

J = 0 

00 DO ll'' I = 1,M 
10 Y { I ) = /, ( I )+ETA+S ( n 
CALL C>-!FCK(Y) 

IF(FAT5.FY.r'''o'''> )CO TO v3 
J = i 

ETA = ET,a /2« 

GO TO 100 

03 IF{J,NF,0)C0 TO - 
EtA = 2 c-KFTA 
GO TO 10+ 

94 E,MN = FtA 
E'„'.v = 2,*FTA 
iro=i 

CALL SEARCH (K, 5) 

.£r'';r.; = n 

e;--ax=e--td 

IC ' = 2 

Call sr,AkcH(x,S) 

ALP;-" =r:'=ID 

kFTsjP.'i 

F i " 

FTC 2''* 

■' 1 1 r 0 j T I r F A R c-j ( X , s :■ 

CO,-; CL/SEPCHi/F+.TNsE:- .6X,F.''!in 

CO.' 'OK'/C-DF/ICD 

cc:. O'- /SERC'-' + /E;-ilNi sE, -AXl 

CC: ON/PCVEL/Z 

CO; ■ .0M/SERCH2/FUNVAL 

CCYHON/SERCHS/F 

EMr!l = FFtI''! 

EMAX1 = F'^'AX 
CCi' .STi=o<.381966 
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CC.,.ST2=..io618i"34 
Dl!'';rr’SlG:\ X(10) »S{inj 
iC DlST = f^'''AX-r'lIM 
DE=GISt 

11 C 1 = F I NS + c 0 M :> T i * D F 
E2 =F'M (. ■'' + C0\>ST2*nF 
Call fd'^'I (x?" .'"IjfI) 
fui=fumval 
Call fmmi {XsS,f2,p2) 

F|J2 = F'IM\/AL 

I'- IF ( ^BS{ F. »0"'1'XAF.5 ( DIST ) ) GO TO 3o 

DF=C0i''i'-.T2--0F 
IF(F1-F2}13,25,14 
14 E;'UW = F1 
El =-F2 
F1 = F2 
Fiu = fu? 

E2=FMINi+COMST2i'-DF 
CALL FiiNl (X,S»E2 ,g2 ) 

Fu 2 = FU'''!VAL 
GO TO 12 
13 EMAX = f^2 
F2=F1 
E2 = E1 
Fu2=FU] 

E1 = FMIF!+C0MST1*DF 
Call Fui'*i(x,s,Fi»ri) 

Ful=FUMVAL 
GO TO 12 

25 IF{ ICD^FO.DGO to 13 
EMI0=F1 

rf.S4X = irp 

DE=FMf.Y-F>',rNi 
GO TO 11 

30 IF(Fl,GT«r2)r-0 TO 32 
FMrO=E3 
Z=F1 
F = F!'1 
F<FTiJRt% 

32 E-;iF = E? 

Z = F2 
F = F'j2 
FCT'JP ; 

T • • 

In;- TO f^L^'l 

Sij-lRC iTT'-tE FU^^Il ( y,S»F,F) 

CCM'.O'.' /i'Ai'jE/s\ 

CO' ""OO/SFRCHA/E I.Xl ,F;ViAXl 
COHLC '/COO:- /I CO 
COi-iLOO/ T'lFCLT ( 30) , SATSFY 
DirFFNSlOi'l X( 10) jS( 10) »Y( iO) 

DO 10 T = l»''i 
in Y( I )=xn ) + -*S( I ) 


IF( TCD«FQ^1)G0 TO 50 
CALL F!.JNT(Y,F! 

i- T"Ri'« ' ■■ 

50 CALL ''■‘-*'"CK( Y ) 

IF.( OAT'^'FYs-JCi. ) GO TO 6o 

F = iOi.>.*< i.-( (E“EMIh 1 )/ (EHAXi-EniNl) )**2y 
kETIJAm: 

6G F = KF + „5 

return 

END 


entry 
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*JOB MEC-24,TlME00ijPAGES0l0»NAME S.iC.CHATURVEDI 

(;^»**KiNe:MaTIC synthesis Of A-BA” fG.< 5-POSlTlON LiHK, r.uinT 


c*** 


SIBJOB 
SIBFTC MAIN 

COMPLEX ONE,F,B ( A) ,C(A),c1(3),ET3{2) ,ET^(2) »EV3{2 ) ,EVA(2 ) ,Z* 


9 


AO 


21 


12 

10 


1Z3»Z4,z5»26 

DIMENSION Tl(8) »RR(8) ,RI (8) ,T2(2) ,V1 (8) ,V2(2) ,G(9) ,U{AI 


DIMENSION X(5),Y(5) 

double PRECISION ZR( lOO ) »ZI ( lOOI ,Sl ( 9 ) 

COMMON/MANE/C,B,T1»T2,v1*V2,X»Y»ONE 

C0MM0N/MAN/g1,G2,g3,G^»H1»h2,H3,hA 

DATA X/10»P»9»6,7,3,5,6,5o3/,Y/11.7,l2e8,llo9,9.7»8»3/,U{l)i 
1U{3) ,U(A)/-12.,-19,,-20.*-19,/ 

PI=A,«-aTAN(1.) 

ONE=( 1. ,0'. ) 

DO 9 1=1,4 
U( n=U( I )*PI/180« 

B( I )=F(l6 ,U( I ) )-ONE 

C( I) =CMPLX(XUI + 1 ) ,Y( I+l) )-CMPLX(X ( 1 ) ,Y( 1) I 
9F{T.E0»1) go to 9 
31(I-1)=C{1)-C( I ) /B{ I )*B( 1) 

continue 

3aLL SYNTH1{G,C1) 

PRINT 40, (G{ I ) ,1 = 1,9) 

F0RMAT(1X,E15<,S) 
call NCRML(G) 

DO 21 1=1,^ 

Sl( n=G(I 5 

Call poly(8,si,2R,zi ) 

DO 12 1=1,8 


RR(n=7R(I) 

Rl (T)=7I ( I ) 

PRINT 10 , (RR(T), 1 = 1 , 8 ), {Ricn»l = l» 8 ) 
FORMAT!///* '"RC ROOTS*, 5 X/ 8 f 15 . 8 //* 


imaginary R00TS*/8F15.{ 


N = n 


DO 11 1 = 1,8 ^ ^ , 

IF(tRI(] '.LT*Oo),OR.(RI(I)»GT.O.)) GO TO 11 


N=N+1 

T1 ( I ) =RP ( I ) 

11 CONTINUE 

DO 13 1 = 1, N 

call SYNTH2( I ,M,T2,T1 ) 

T 1 ( n= 2 .*ATAN(Tl ( I) ) 

IFfMcEQoO) GO TO 13 
DO 13 J=1,M 

ET3(J)=(G1*F(1o,T1(I) )+G2*F( 1. ,T2( J) )+G4)/63 
Et 4( J) = (Hl*F ( ‘ . sTl ( n ) +H2*F(*. ,T2 ( J ) )+h4)/h3 
call SYNTH4( I ,J,Z2,Z5,Z6) 

C(l}=Z 2 *{F{ lo ,T1 ( I ) )-ONE) 

C(2)=Z 2 *(F(1. ,T2(J) )“ONE) 
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C(3)=Z 2 *(ET3{ J)-ONE) 

C(4)=Z 2 *(ET4( J)-ONE) 

DO 14,K=1,4 
B{K)=^B(K) 

IF(K.EQ*1) GO TO 14 

Cl (tC-l )=C( l)-C(K) /6(K)*B( 1) 

14 continue 

3aLL SYNTH1(G,C1) 

Print 40» (G( im) r,iM=i,9) 

Call nofml(G) 
do 22 IM=1,9 
22 S1{IM)=G(IM) 

Call. poly(8»si,2R,zi ) 

DO 16 L=1,8 
RR(L)=7P(L) 

16 Rl(L)=7nL) 

print 17,{RR<L) ,L = 1,8) ,(RnL),L=l,8) 

17 FORMAT!///* REAL R00TS*»5X/8f 16,8//* IMAGINARY ROOTS^/FFlS.f ) 
Nl=o 

DO 18 L=1.8 

IF { (RI (L' oLT.Oe. } .OR. ( RK L) ,GToO. ) ) GO TO 18 

Nl=Nl+l 

V1{L)=RR(L) 

18 continue 

DO 13 L=1,N1 

Call synth2(L,mi,v2,vi ) 

VI {L)=2.*ATAN(V1 (L) ) 

IF(M1,EQ«0) GO TO 13 

DO 13 !<;=1,M1 \ 

EV3(K) = (G1*F( 1. ♦VKL) ) +G2*F ( 1 . , V2 { K ) )+g4)/g3 
EV4(K)=(H1*F(1. »V1(L) )+H2*FC 1.»V2(K) )+H4)/h3 
Call synth5{l,k, z4,z3,zi,z2) 
print 30,Z1,Z2,Z3,Z4,Z5,Z6 
30 F0RMAT(//11X.6(F15.8,5X) ) 

13 CONTINUE 
STOP 

end 

sibftc synthi 

subroutine synthi {6»v1) 

dimension G(9),T!(8),T2{2),V*{a),v2<2),y(;.,,Y{-'} 

COMMON /M an /G i *62 , g3 , GA , H 1 , , K-> 9 H4 

COMMON/MANE/C »d *71 »T2 , Vl ,V2 , X » Y ,0NE 
>.wKM0N/S0B4/A1,A2 , A3 , a4 , a5 , a6 , B1 > B2 , B3 , B-. 

COMPLEX P^QjR.S »P1»Q1 ,R1 ,S1»0NE»Z.- 

COi-iPLEX B(4) »C{4) 5Yl(3)»Gl?G25G3»G49G*'»Hi»H.‘>H«-»H“‘*n^ 
Gl=B(3i*B(25*'ONE-Yl{l)/Yl(2) ) 

G2=-B(3)*B( 1 ) 

G3=e( 1 )*B( 2)*Y1 (1) /Y1 ( 2 ) 

G4““* ( G1 + G2+g3 ) 

G5=Gl*rONJG ( Gl ) +0-*Cv»<wGlG2}+r;4«-rO!^J5{G4)-G3*CONJG(G3 ) 
Hl=B(2)*B(45*(ONE--Yiil J /Y’ (3) ) 

H2®-B ( 1 3 Y i 

1 5*B( ?)*Yl (1 )/Yl { 3) 
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H4=-(H1+H2+H3) 

H5=Hl*rONJG{ Hi )+H2*C0NJG(H2 )+H^*C0NJG(H^)--H3*C0NJG{H3 ) 
P=61*C0NJG(G2 ) 

Q=G1*CCNJG(64) 

S=G2*C0NJG(G4i 
R = G5 

P1=H1*C0NJG( H2 ) 

Q1=H1*C0NJG< HA) 

Sl=H2*c0fMJG{H4) 

R1=H5 

A1=REAL(P-Q+R-S) 

B1^-2,«AIMAG(P+Q) 

C1=REAL(Q-P+R-S) 

A3=REAL( R+S-P-Q) 

B3=2.*AIMAG(P-0) 

C3=REAL(P+Q+R+S} 

A5=2o*AIMAG(P-S) 

B5=4,*REAL (P ) 

C5=-2*^fAlMAG(P+S) 

A2=REAL( Pl-Ql+Rl-Sl ) 

B2=-2.*AIMAG(Pi+Ql5 

C2=REAL(Q1-P1+R1-S1) 

A4=REAL(R1+51-P1-Q1 ) 

B4=2.*aIMAG(P1-Q1) 

C4=REAL( Pl+Qi+Rl+Sl ) 

A6=2.*aIMA6(P1“S1 ) 

B6=4,*REAL{P1 ) 

C6=-2**AIMAG(P1+S1} 

D1=A4*a5-A6*A3 

02 =A4*p5+a5*B4-a3*86-A6*B3 

D3=a4*c5+a5*C4+84*B5-a3*C6-A6*C3-S3*B6 

D4 = B4^«-c5+B5*C4-B3*C6-B6*C3 

D5=C4*c5-C3*C6 

E1=A1*a6~A2*A5 

E2=a1*b6+A6*B1-a2*B5-a5*b2 

E3=a1*C^+a6*C1+b1*B6-a2*c5-a3*C2-s2*b5 

E4=Bl*r6+B6#Cl-B2*c3-35*C2 

E5=C1*C6“C2*C5 

F1=A1*a4-a2»A3 

F2=A1*B4+A4*81-A3*B2-A2*B3 

F3=A 1*g4+a4'«-C1+B1*B4-a3*-c2-a2*C3-B3*82 

F4=R4*c1+B1*C4-B^^C2-B2*C3 

F5=C1*g4-C2*C3 

6( l)=D5*E5-F5sf--2 

G( 2) =D4«-E5+D5*E4-2,*F4*F5 

G ( 3) =D3*E5+D4*E4+D5*E3-2,^cP3^tF5-F4**2 

G( 4) =D2*E5-4-D3^!-E44-D4*E3+D5i^E2-2o*F2*F5-2,*F3*F4 

G ( 5 ) =D1*E5+D2*E4+L')3a-F3+D4*E2+D3*E 1-2 <,*f1*F5-2,-k-F2*F4— f 3**2 
G{ 6)=Dl’^E4+D2*E3+n3*E2+D4*El-2,*Fl*F4-2,#F2*F3 

G< 7)=Dl*E3+D3*‘^l+n2^^E2-2.*Fl*F3-F2«-*2 

G{ 8)=D2*E1+D1*E2-2,*F1*F2 
Gf 9)=D1^E1-F1»*2 
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RETURN 

End 

SIBFTC SYNTH2 

SUPROUTINF SYNTh2(I ,M,T2,T) 

LOf'’’.'’iO!'M/SUFs2/ a1 9 A 2 sA3 >A^ > A6 jBI »92 

DInENSiGN T2(2 ) 9XR(2 ) »xn 2) 9YR(2 ) jYI (2) ,X(2 ) »Y(2 ) »RR( 3 ) »RI (3) ,T(8 ) 
RR { 1 ) =FIJf.;c ( A3 , 3 , c3 , T , T ) 

RR ( 2 ) =FUNC ( A3 , 35 ,r3 , f , t ) 

RR ( 3 ) =FUNC ( a i 9 U , Cl 9 1 5 T ) 

R T { 1 ) =PUMC ( A '= 9 BA , CA , I T ) 

Rl {2)=FUNC(AB*B69C6,T) 

Rl(35=~i!XC(A.?.Y2,r2,T) 

Call quad (RR , 'R ,x i) 

Call GijAD(^^i ,yR,y! ) 

M = 0 

DO 14 J=l,? 

IFnxi(J)»LT.Oo)«OR,(XI(J),GT.C.)) RETURN 
X{J)=XR( J) 

15 IF( (YT ( J) .LT.n. ) «OR, (YI ( J)«GT*0. ) ) KFTURN 
Y( J)=YR( J) 

14 continue 

DO 16 J=l,2 

IF{<X{J).GT,Y(J) ).OR. (X( J)*LT.Y(J} ) ) GO TO 15 
T2( J)=2.*ATAN( XCJ)) 

M=M +1 

16 CONTINUE 
IF(f'’.EQ.2) RETURN 
P=X( 1 ) 

.Q=X(2) 

DO 17 J=l ,2 

IF ( (X ! J) .GT.Yt J) ) .OR. (X( J) ,LT.Y{ J) ) ) GOTO 17 
Y=M +1 

T2 ( J)=2.-^ATAM(X (J)) 

17 continue 
RETURN 

end 

SIBFTC SYNTH4 

Subroutine synthac i ,J9z2 ,z5,z6 ) 

C C ^ XM c ' 7 '• A N F / r 9 B » T 1 9 T 2 , V 1 , V 2 9 X » Y 9 0 N E 
D I mens T ON X ( 3 ) , Y ( 5 ), , T 1 ( 8 ) 9 T2 [ 2 ) , V 1 ( 6 ) 9 V2 ( 2 } 

CCFPLrx Z6,Z2,25 ,c11 9A(4) ,8(4) ,C{3) ,F90NF9A11 
A(i)=P(i.,Tl{I5)-ONE 
A (2!=:7 1 • ,T2 f J) )-ONE 
A 1 1 = A { 1 } -A { 2 5 /B ( 2 ) ( 1 ) 

CJ 1=C(1 }-C(2)/B{2) 

25 =C11/A{1) 

Z2 =(r(l)-3{l)*Z 5 )/A{ 1) 

Z 6 =F{X(1) 9Y{1 ) )-22-Z5 

RETURN 

End 

TIBFTC SYNTH5 

subroutine S • i T H 5 { L , , Z 4 , Z 3 , 2 1 9 Z 2 ) 

COKFON/Y aNp/C ,R 9 T 1 , T2 , Vi ,^2 , X , Y ,ONF 


dimension Xt5)5Y(5),Tl(8)5Vl{8),T2(2)?V2{2i 
COMPLEX B(4) sC{3) ,Zl , Z 2 , Z3 , Z4 , Cl 1 , F , ONE » A ( 4 ) , A 11 


A ( 1) =F{ 1= ,V1 (L ) )-ONE ■ 
A(2)=F(1o »V2{K) }-0NE 
A11=A(1)-A(2)/B(2)*B(1) 
Cll= C{ 1 )-C{ 2 ! /8 ( 2 ) 
Z3a=Cll/A ( 1 ) 

2A=(C{1)‘^B(1S*Z5WA<1) 

Z1=Z2+73“ZA 

t^c TURN 

tixD 

ffTopfr friiNc 


FUNCT I ON FUNC ' ,0 , S * T ) 

COMHON/MANE/C ,3 » T I , T2 , Vi ,V2 , X ,Y,ONr 

Dimension tks) ,t 2(2) »vic8) ,v 2(25 ,x( 5 ) 

COMPLEX C( 3 ) ,B( 4) ,0NE»Z2 

FUNC=P*T { I )*-^2+Q*T( I >4-S 

Return 

END 

$I6FTC QUAD 

Subroutine quad(r,x,y) 

Dimension R(3) ,x(2 ) ,y(25 
X(1)=-R(2)/(2.*R(1) ) 

X<2)=X{1) 

RADCL = R ( 2 ) **2-4.*R ( 1 ) -SJ-R ( 3 ) 

IF(RADCL) 2o»21,22 

20 Y( 1)=SQPT(-RADCL ) /(2,*R{ 1) ) 

Y(2)=-Y{1) 

RETURN 

21 Y(1)=0 
Y(2)=0 

22 X(1)=X{1'+SQRT(RADCL)/(2.*R( 1) ) 
X(2)=X(1)-SQRT(RADCL)/(2,*R( 1) ) 

Y{1)=. 

Y(2)=:r> 

return 

END 

SIBFTC F 


»Y(5) ,T(9) 


COMPLEX FUNCTION F(X»Y) 
F=cKPLx{X, 0 . )*CMPLX(C0S{Y) ^SlNtY) ) 
RETURN 

end 

SIBFTC NORML 

Subroutine normlig’ 

Dimension g(9) 

X=ABS{Gf 1 ) ) 

DO 10^ 1=2,9 

if{x.gt.abs(G( n ) ) x=G(i) 

100 continue 

DO lol 1=1,9 

101 G(n=G{I)/X 
Return 


END 

■ETCH POLY MEG134 

c POLY IS Library subroutine at i,ieTe kanpur computer centre 
c and has been called from the disk 
c subroutine poly (NM,CC,RR»Rn 

c MAIN LlriE FOR COKPUT I NS- ROOTS , R , OF THE POLYNOMIAL, 

c CC<l)+rr (2)*X+CC(3)*(X*^2)+ , +CC ( N+1 ) * ( > 

C ROOTS APE CCMRUTED FROM THE COFFICIENTS VIA THE MULLER METHOD 
C N E-WTO Nc METHOD. 

C RR(NN) ARE THE NN REAL ROOTS OF THE POLYNOMIAL* 

C RI(NN) >RE THE NN IMAGINARY ROOTS OF THE POLYNOMIAL, 

DOUBLE PRECISION tP1,RR( 100) sRlI iOO) >CC{9) »CCC{l0l) ♦SR(S) 

>ONE 

rrsTRY 


j 111 . ■■ ■ 

iGBlifr- ■■■ '■ -’AT/aY, 


i 






REFINED BY 

POLCOl 



